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Abstract 

The adsorption and dissociation of O2 molecules at the Be(OOOl) surface is studied by using 
density-functional theory within the generalized gradient approximation and a supercell approach. 
The physi- and chemisorbed molecular precursor states are identified to be along the parallel and 
vertical channels, respectively. It is shown that the HH-Z (see the text for definition) channel is the 
most stable channel for the molecular chemisorption of O2. The electronic and magnetic properties 
of this most stable chemisorbed molecular state are studied, which shows that the electrons transfer 
forth and back between the spin-resolved antibonding vr* molecular orbitals and the surface Be sp 
states. A distinct covalent weight in the molecule-metal bond is also shown. The dissociation of 
O2 is determined by calculating the adiabatic potential energy surfaces, wherein the T-Y channel 
is found to be the most stable and favorable for the dissociative adsorption of O2. Remarkably, we 
predict that unlike the other simple sp metal surfaces such as Al(lll) and Mg(OOOl), the adiabatic 
dissociation process of O2 at Be(OOOl) is an activated type with a sizeable energy barrier. 

PACS numbers: 68.43.Bc, 82.20.Kh, 82.45.Jn, 34.80.Ht 
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I. INTRODUCTION 



It is of great scientific importance to understand the beliaviors of diatomic molecules at 
solid surfaces, including their adsorption and dissociation, and corresponding energy barriers 
during the bond breaking and bond formation at the surfaces l2|. Of all prototypes, the 
interaction of O2 molecules with metal surfaces has gained lots of interest for many techno- 
logically relevant processes such as heterogeneous catalysis and corrosion [3]. Theoretically, 
ab initio modeling has been successfully used over a wide range to study the adsorption and 
dissociation of O2 molecules at transition metal surfaces. By calculating the adiabatic poten- 
tial energy surface (PES), it has been found that O2 molecules will spontaneously dissociate 
while adsorbing at reactive transition metal surfaces like iron (te 

l9| and Nickel (Ni) [9^, the 



■on (F( 
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metals like gold (Au) [5|], silver (Ag) , platinum (Pt) 
adsorption of O2 turns out to depend on the ambient temperature. Remarkably, in all above 
transition metal systems, the concept of adiabatic PES works very well in explaining and 
predicting a large amount of physical/chemical phenomena during dissociation process of 
O2. 

When the attention is focused on the simple sp metals such as Al(lll), however, an un- 
comfortable gap opens between ab initio prediction and experimental observation. The most 
notable is the long-term enigma of low initial sticking probability of thermal O2 molecules 
at Al(lll), which has been measured by many independent experiments Q, [ill tiut cannot 
)e reproduced by adiabatic state-of-the-art density functional theory (DFT) calculations 



12, [l^- The central problem is that the adiabatic DFT calculations were unable to 
find any sizeable barriers on the adiabatic PES, whose presence, however, is essential for 
explanation of the experimental finding. This has led to speculations that nonadiabatic 
effects rnay play an important role in the oxygen dissociation process at the Al(lll) surface 
[13I . 15, 16, 17, 18, Q, Recently, a semiquantitative agreement with the experimen- 



tal data is achieved by nonadiabatically confining the trajectories of the approaching O2 
molecules to the spin-triplet PES 2l|] . Until now, however, it still remains unclear how and 
in what chemical circumstance this spin-selection rule is reasonable, and different opinions 
exist in literature 2^- For oxygen dissociation at another simple metal, Mg(OOOl), first- 



principles DFT calculation also shows the lack of a barrier on the adiabatic PES [23] , which 
is in strong disagreement with experimental observation 24, l25| of a low sticking coefficient 
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of O2 at Mg(OOOl). This discrepancy was ascribed, in a similar manner with that in the 
02/Al(lll) system, to nonadiabaticity in the dissociation process of oxygen molecules 23|. 

Therefore, it becomes clear that the theoretical study of oxygen dissociation at the simple 
sp metals is still far from its maturity. In particular, considering the fact that up to now most 
of the theoretical results and conclusions are only based on the 02/Al(lll) and O2/Mg(0001) 
prototypes, workers are thus confronted with an important question: Does it retain necessary 
or valid for all simple sp metals (at least, for metals with similar elemental valence electrons 
to Al or Mg) to take into account such nonadiabaticity as spin-selection rule in obtaining an 
activated PES? Motivated by this question, in this paper we have carried out a systematic 
ab initio investigation of the adsorption and dissociation of oxygen molecules at Be(OOOl) 
surface. Beryllium has the same crystal structure and valence electrons as magnesium does. 
Subsequently, it is attempting for one to derive that the behaviors of oxygen dissociation 
at Be(OOOl) are the same as or similar to that at Mg(OOOl). Our results, however, show 
that this is not true. The most distinct is that in the present O2/Be(0001) system, our 
calculated adiabatic PES displays sizeable energy barriers along the dissociation paths. This 
partially but definitely answers the above question. That is, in the DFT calculations of the 
simple-metal systems, the inclusion of nonadiabatic effects is not always indispensable for 
the presence of molecular dissociation barrier. 

Besides this basic point of interest, our present study is also motivated by the fact that Be 
has vast technological applications due to its high melting point and low weight. During these 
applications, surface oxidation as the main kind of corrosion always needs to be prevented. 
Thus a systematic study on the adsorption and dissociation of O2 molecules at Be surfaces, 
i.e., the initial stage during the surface oxidation process, should be done. Moreover, Be 
is also a getter in experimental nuclear fusion reactors to adsorb residual gases such as O2 
and H2O in the plasma vessel, improving the plasma cleanliness 26|, which also makes it 
highly meaningful to study the adsorption and dissociation of O2 at Be surfaces. At present 
there are only very few experimental data on the oxidation of Be, reporting that the surface 
oxidation begins at an oxidation nucleation center, followed by a spreading growth from the 



center 
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291] , while ab initio studies are entirely lacking. This as a whole encourages 



us to theoretically report a systematic investigation on the O2/Be(0001) system. The rest of 
this paper is organized as follows: In Sec. II the computational method is briefly described. 
In Sec. Ill results for the free oxygen molecule, the bulk Be and clean Be(OOOl) surface are 
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given. In Sec. IV and Sec. V we present our results for O2 adsorption and dissociation 
at Be(OOOl) surface, respectively. Finally, we close our paper with a summary of our main 
results. 



II. COMPUTATIONAL METHOD 



The density-functional theory (DFT) total energy calculations were carried out using 
the Vienna ab initio simulation packa ge [sO ] with the projector-augmented-wave (PAW) 



31| and plane waves [3^. The so-called "repeated slab" geometries were 



pseudopotentials 



employed 



33l |. This scheme consists in the construction of a unit cell of an arbitrarily fixed 
number of atomic layers identical to that of the bulk in the plane of the surface (defining the 
dimensional cell), but symmetrically terminated by an arbitrarily fixed number of empty 
layers (the "wacunm") along the direction perpendicular to the surface. In the present 
study, the clean p{2 x 2)-Be(0001) surface was modeled by periodic slabs consisting of nine 
Be layers separated by a vacuum of 20 A, which was found to be sufficiently convergent. 
Oxygen molecules were symmetrically introduced on both sides of the slab. During our 
calculations, the positions of the outmost three Be layers as well as the O2 molecules were 
allowed to relax until the forces on the ions were less than 0.02 eV/A, while the central 
three layers of the slab were fixed in their calculated bulk positions. The plane-wave energy 
cutoff was set 400 eV. After a careful convergence analysis, we used allxllxl fc-point 
grid for the p{2 x 2) cell with Monkhorst-Pack scheme 3J|. Furthermore, the generalized 
gradient approximation (GGA) of Perdew et al. 35|] for the exchange-correlation potential 
was employed since the GGA results has been previously validated for the bulk Be 36|. A 
Fermi broadening 37| of 0.1 eV was chosen to smear the occupation of the bands around 
Ep by a finite-T Fermi function and extrapolating to T=0 K. 



III. BULK BE, CLEAN BE(OOOl) SURFACE, AND FREE OXYGEN MOLECULE 



First, the total energy of the bulk hep Be was calculated to obtain the bulk lattice 
parameters. The calculated lattice constants are 2.26 and 3.56 A respectively for a and c, 
according well with experimental values 38| of 2.285 and 3.585 A. The band structure and 
orbital-resolved density of states (DOS) of bulk Be is shown in Fig. 1. Typically, the DOS 
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FIG. 1: (Color online). Band structure (left panel) and orbital-resolved DOS (right panel) of the 
bulk hep Be. The Fermi energy is set at zero. The inset shows the Brillouin zone. 

of bulk Be has a relatively small value at the Fermi energy [Ep), because there exists a 
Dirac point at the Fermi energy in the band structure along the F— M direction, and a wide 
band gap along the F— A direction. This anisotropic property of Be differs from that of 
other alkali-earth metals, whose band structures around the Fermi energy are always nearly 
free-electrons like. Although the electronic configuration of elemental Be is ls^2s^, one can 
see from Fig. 1 that the 2p states play an important role in the DOS behaviors below Ep, 
indicating a strong hybridization between s and p electronic states in bulk Be. 

The band structure and total DOS for the clean Be(OOOl) film are shown in Fig. 2. 
Compared to the bulk case (Fig. 1), one can see from Fig. 2 that there exists no longer 
Dirac points or band gaps along any direction. Moreover, the DOS at Ep is prominently 



enhanced. This result accords well with previous studies [3; 



Similar results have also 



been observed for W (110) and Mo (110) films |4ll]- From the inset in Fig. 2, we can see 
that the surface electronic states around Ep mainly accumulate within the two topmost 
Be layers. Further detailed wavefunction analysis shows that these states are mainly Be 
2p states. Our clean-surface calculation shows that due to this pronounced surface charge 
redistribution, the two outmost Be(OOOl) layers relax significantly from the bulk values. The 
first-second interlayer contraction is 3.8% and the second-third interlayer expansion is nearly 



1.2%, which is in agreement with recent first-principles calculations 421] and comparable with 
experimental measurements 

The total energies of the isolated O atom and free O2 molecule are calculated in an 
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FIG. 2: (Color online). Band structure (left panel) and total DOS (right panel) of the clean 
p{2 X 2) Be(OOOl) surface with nine atomic layers included in the supercell. The Fermi energy is 
set at zero. The inset shows the vertival distribution of the charge density close to the surface 
(Nos. 1-4 enumerate the atomic layers from the surface). 
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FIG. 3: (Color online). The calculated spin-polarized density of states for the molecular orbits 
(MO's) of the free O2 molecule. The Fermi energy level is set at zero. 

orthorhombic cell of scale 13 x 11 x 17 with a (3 x 3 x 3) /c-point mesh for the Brillouin 
zone sampling. The spin-polarization correction has been included. The binding energy of 
O2 is calculated to be 1/2E°2=2.89 eV per atom and the 0-0 bond length is about 1.235 
A. These results are typical for well-converged DFT-GGA calculations. Compared to the 



experimental 



44 1 values of 2.56 eV and 1.21 A for O binding energy and bonding length, the 



usual DFT-GGA result always introduces an overestimation, which reflects the theoretical 
deficiency for describing the local orbitals of the oxygen. 

In order to show the general charge-transfer or redistribution effects on O2 molecular 
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orbits (MO's) after adsorption, here we first calculate and plot in Fig. 3 the MO-resolved 
local density of states of the free O2 molecule. Typically, the bonding a MO's are lower 
in energy than the bonding vr MO's, for both spins. This ordering of the MO's has im- 
portant consequences for the molecular bond to the surface. Also it reveals in Fig. 3 that 
spin splittings are sizeable (~2 eV) for both bonding and antibonding MO's. This reflects 
Hund's spin rule, which describes ground state as spin-polarized for O2 (S'=l). The highest 
occupied MO (HOMO) and the lowest unoccupied MO (LUMO) are the spin-up and spin- 
down antibonding vr* MO's, respectively. Our calculated picture of t^ MO's for a free O2 
molecule is in good agreement with the previous theoretical reports 
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IV. THE ADSORPTION OF O2 MOLECULE 

There are four high-symmetry adsorption sites on the Be(OOOl) surface, including top (T), 
hep hollow (HH), fee hollow (FH), and bridge (B) sites depicted in Fig. 4(a). In this study, 
we construct twelve initial structures with high symmetries by orienting the O2 molecule 
at the four high-symmetry sites respectively along the X (i.e., [1120]), Y (i.e., [1100]) and 
Z (i.e., [0001]) directions. We also construct several low-symmetry initial structures by 
rotating the O2 molecule in the XY , YZ , and XZ planes with small angles. In all the initial 
configurations, the heights of the O2 molecules are set at /io=4 A, see Fig. 3(b). It is found 
within our expectation that after geometry optimization, these low-symmetry structures will 
either relax into the high-symmetry ones or be less stable with a lower adsorption energy 
than the high-symmetry ones. This is similar to what has been observed in studying the 
adsorption of O2 molecules at Pb(lll) [47]. After geometry optimization, it is found that all 
the parallel adsorption states with the O2 molecule lying down on the substrate surface are 
stable. For the vertical entrances with the O2 molecule oriented perpendicular to the surface, 
one exception is the B-Z entrance (namely, the O2 molecule at the bridge site with the 0-0 
bond along the Z direction), which turns to evolve into the HH-Z adsorption configuration 
after geometry optimization. 

One central quantity tailored for the present study is the average adsorption energy of 
the adsorbed oxygen species defined as 

1 1 

Ead = — ^r^[-£^O2/Be{0001) " -E'Be(OOOl) — Nq X --E02]) (l) 
iVo ^ 



7 




FIG. 4: (Color online) (a) The p(2 x 2) surface cell of Be(OOOl) and four on-surface adsorption 
sites. Here only the outmost two layers of the surface are shown, (b) The sketch map showing that 
the molecule (with vertical or parallel orientation) is initially away from the surface with a hight 
/iQ. Note that although the atomic radius of O is much larger than that of Be, for convenience 
in depicting the adsorption positions, here and in the following we use larger circles to denote Be 
atoms. 

where A^'o is the total number of O atoms presented in the supercell, -E'o2/Be{oooi)) -E'bc{oooi)5 
and E02 are the total energies of the slab containing oxygen, of the corresponding clean 
Be(OOOl) slab, and of a free O2 molecule respectively. According to this definition, a positive 
value of Ead indicates that the adsorption is exothermic (stable) with respect to a free O2 
molecule and a negative value indicates endothermic (unstable) reaction. 

Starting from the above mentioned entrances and after geometry optimization, the ob- 
tained molecular adsorption energy (Ead), molecular magnetic moment (MM), work function 
($), adsorption height {h), and 0-0 bond length {d) are listed in Table I. From Table I the 
following molecular adsorption features are revealed: (i) For the O2 molecule in parallel with 
Be(OOOl) surface, the relaxed molecular structures vary very little compared to the initial 
molecular structures, which suggests the molecular adsorption to be the physisorption in 
these parallel channels. In these cases, the relaxed adsorption height is around 3.9 A. The 
0-0 bond length is 1.236 A, displaying a negligible expansion when compared to that of 
a free O2 molecule (1.235 A). The molecular MM almost saturates at its free value of 2.0 
/ifi, also suggesting the unaffected molecular orbitals in these physisorbed structures. In 
addition, the work function of the adsorbed Be(OOOl) surface is almost identical to that 
of the clean Be(OOOl) surface, implying negligible charge transfer between the parallel O2 
adsorbate and the surface Be atoms. Finally, the calculated molecular adsorption energy is 
only ~23 meV, which is so small that a little thermal fluctuation may result in desorption 
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of the O2 adsorbate with parallel 0-0 bond from the Be(OOOl) surface. This result is con- 
sistent with the recent experimental observation that the reactivity of oxygen atoms with 
Be(OOOl) surface is greater than that of oxygen molecules with Be(OOOl) surface 29|. On the 
other hand, one may question whether, when overcoming the energy barrier by artificially 
setting the initial metal-molecule distance Hq to be lower than the above obtained physisorp- 
tion height (~3.9 A), new kinds of molecular adsorption state with parallel 0-0 bond will 
occur or not. For this, a large amount of calculations with more lower initial adsorption 
heights have been carried out and no other parallel molecular physi- or chemisorption states 
have been found, for details see discussion in the next section; (ii) For the O2 adsorbate 
with the 0-0 bond oriented perpendicular to Be(OOOl) surface, as shown in Table I, the 
relaxed molecular structures display a signature of chemisorption. Among them the HH-Z 
configuration is most stable with the largest molecular adsorption energy of ~0.5 eV (per 
atom). In this channel the adsorption height is decreased to be 1.58 A, and prominently, 
the molecular MM is largely decreased to be 0.8 /is, which suggests a pronounced charge 
redistribution among the MO's via interaction with the Be(OOOl) surface. The 0-0 bond 
in this most stable chemisorption channel is 1.471 A, indicating a large expansion from 
the free O2 molecule and a fundamental weakening of the molecular bonding. In addition, 
the work function change is also prominent, implying an observable charge redistribution 
between the chemisorbed O2 molecule and the Be(OOOl) surface. Meanwhile, during the 
chemisorption of the O2 molecule, the Be(OOOl) surface is also influenced. Specifically, the 
three Be atoms around the adsorbed O2 molecule of the HH-Z entrance are pulled out by 
about 0.3 A. Recalling that the bulk BeO has an unusual wurtzite structure, it is within 
one's expectation that O2 chemisorption in the HH-Z channel (instead of the FH-Z channel) 
is most energetically favorable. Our additional calculations of the atomic oxygen adsorption 
at Be(OOOl) also show that the hep hollow site is most stable for atomic oxygen adsorption. 

To clarify the bonding interaction between the chemisorbed O2 molecule (in the HH-Z 
channel) and the Be(OOOl) surface, we calculate and plot in Fig. 5(a) the electron-density 
difference Ap(r), which is obtained by subtracting the electron densities of noninteracting 
component systems, PBe(oooi) (i") "I" P02('")) from the density po2/Be(oooi) 

(r) of the O2/Be(0001) 

system, while retaining the atomic positions of the component systems at the same location 
as in O2/Be(0001). A positive Ap(r) obviously represents charge accumulation, whereas a 
negative Ap(r) represents charge depletion. One can see from Fig. 5(a) (where the contour 
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TABLE I: The calculated adsorption energy per atom (Ead), molecular magnetic moment (MM), 
work function (<I>), adsorption height (h), and 0-0 bond length (d) for physi- and chemisorptions 
along different channels. 



Channel Ead (meV) MM (/xs) h (A) d (A) $ (eV) 
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spacing is 0.02e/A"^) that the charge redistribution mainly occurs at the surface and involves 
the chemisorbed O2 molecule and the two topmost Be (0001) layers. It is apparent that 
upon molecular chemisorption, the occupation of the surface and subsurface Be sp states 
is decreased, while the occupation of O2 tt* antibonding MO's is increased, which suggests 
an electron transfer from the former to the latter. Meanwhile, it reveals in Fig. 5(a) that 
the occupation of other bonding (a, tt) and antibonding MO's (a*) is also decreased. The 
target states for these transfered MO charges are of course again the O2 tt* antibonding 
MO's, as shown in Fig. 5(a). A well-known harpooning mechanism is responsible for this 
inter-MO charge transfer. Prominently, the electron redistribution of the O2 tt* antibonding 
MO's is largely asymmetric with respect to the two O atoms. The reason is that the lower 
O atom in the O2 molecule is in a complex bonding state. Besides interacting with the 
upper O atom to form a weakly-bound molecule, this lower O atom also strongly interacts 
with the surface and subsurface Be atoms via a mixed ionic/covalent bonding mechanism. 
In particular, it is the strong covalcncy in the Be-0 bond that heavily distorts the charge 
distribution of the O2 tt* antibonding MO's. To our knowledge, this strong covalency in the 
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FIG. 5: (Color online) Contour plots of the (a) difference electron density and (b) the spin density 
for the relaxed O2/Be(0001) slab with O2 molecule chemisorbed along the HH-Z entrance. Solid 
and dotted lines denote accumulated and depleted densities, respectively. 

bond between the metal surface and oxygen molecule has not been theoretically reported in 
previous studies. On the whole, the charge transfer to the vr* antibonding MO's, as well as 
the strong covalency in the Be-0 bond, weaken the molecular bond and spin polarization. 
For spin density of the chemisorbed O2 molecule in the HH-Z channel, see Fig. 5(b). 

In order to gain more insights into the precise nature of the chemisorbed molecular state 
in the O2/Be(0001) system, the orbital-resolved site-projected densities of states (PDOS) for 
the O2 molecule (in the most stable HH-Z entrance) and the topmost Be layer are plotted 
in Fig. 6(a) and 6(b), respectively. By comparison with the case of a free O2 molecule, 
one can see that the MO properties of the chemisorbed O2 in the HH-Z entrance undergo 
the following fundamental changes: (i) The spin-split PDOS signature (peaks) for the two 
energy-lowest bonding and antibonding a MO's for both spins in Fig. 3 changes to vanish 
in Fig. 6(a), indicating their breaking in interacting with the Be(OOOl) surface. This is in 
accord with the electron-density difference result shown in Fig. 5(a) that charges mostly 
flow out of these two MO's; (ii) The PDOS peaks for the two nearly-degenerate bonding 
(Tp and TTp MO's around E=—7.5 eV are broadened upon chemisorption to merge together, 
with the amplitude becoming much weaker than that in free O2. Also, the spin splitting of 
these two bonding MO's vanishes. These features consistently reveal that although these 
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FIG. 6: (Color online) The orbital-resolved PDOS for the chemisorbed O2 molecule along the HH-Z 
channel and (b) for the topmost Be layer. The Fermi energy is set at zero. Note that for O2 the 
PDOS is plotted in a spin-split form. 

two bonding MO's do not lose their nature as the MO's, they are largely reshaped due to 
their hybridization with the sp orbitals of the surface and subsurface Be atoms. This metal- 
molecule hybridization is so strong that there develops, as shown in Fig. 6, a sharp sp-hybrid 
peak around E=—7.5 eV in the PDOS of the topmost Be layer; (iii) The spin-up antibonding 
71* MO shifts up in energy towards Ep and is partially depopulated by an observable amount 
of charges. As a compensation, the spin-down antibonding it* MO shifts down towards Ep 
and becomes partially occupied. Obviously, the spin rule prohibits direct charge transfer 
from the spin-up to spin-down antibonding vr* MO. Therefore, the driving or intermediate 
factor for this charge decrease in the spin-up and increase in the spin-down it* MO's is due 
to the metallic Be(OOOl) surface, which accepts electrons from the spin-up vr* MO and then 
donates electrons to the spin-down tt* MO. From this aspect, the present result is consistent 
with the well-known harpooning mechanism, which describes adsorbate-metal interaction by 
a simultaneous transfer of electrons from the adsorbing molecule into the unoccupied metal 
states (direct bonding) and a back donation of electrons from occupied metal states into 
the antibonding adsorbate orbitals. In the present case, this exchange is achieved through 
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FIG. 7: (Color online) Color filled contour plot of the potential-energy surface for O2 dissociation 
at Be(OOOl) along the T-Y channel, as functions of the O2 bond length d and distance h (from the 
surface). The contour spacing is 0.1 eV. The dashed line indicates the reaction pathway. 

interaction with the sp orbitals of the Be(OOOl) surface; (iv) Even more interestingly, to 
compensate the breaking of the two energy-lowest bonding and antibonding a MO's, there 
develops a new hybrid peak at £'=—22.7 eV in the PDOS of the system. From Fig. 6 one 
can see that this peak is characterized by a strong hybridization of O 2s and Be sp states. 

V. THE DISSOCIATION OF O2 MOLECULES 

In order to deepen our understanding of the initial stage of oxidation at Be (0001) surface 
from a theoretical point of view, we also calculate the potential energy surfaces (PES's) 
for O2 molecules on the Be (0001) surface, which depicts the adiabatic dissociation path 
with the lowest energy barrier. Since it has been shown that the PES given by the density 
functional theory gives the result that can be compared with experiments, we expect that 
our calculated PES provides the qualitative feature of the molecular dissociation process in 
the present specific O2/Be(0001) system. Figure 7 shows the obtained PES as a function of 
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the distance h from the surface and the bond length d of the O2 molecule along the T-Y 
channel, which, after a large amount of calculations, turns out to be the most favorable 
dissociative adsorption channel with largest adsorption energy (4.14 eV/atom) and lowest 
activation barrier (0.23 eV). 

One can clearly see from Fig. 7 the process from initial molecular physisorption with a 
height of 3.9 A and bond length of 1.23 A to the final dissociative adsorption. The most 
distinct feature in Fig. 7 is that the calculated adiabatic PES exhibits a sizeable energy 
barrier toward dissociative adsorption. The transition state for this barrier is located at a 
distance of 2.2 A from the Be surface, where the 0-0 bond length is elongated to 1.30 
A. Correspondingly, if we initially put O2 in the T-Y channel at a height lower than 2.2 
A, then do geometry optimization, the O2 molecule will spontaneously dissociate with the 
two O atoms moving to the hep hollow and fee hollow sites, respectively. Considering that 
the previous ab initio calculations have failed to obtain a barrier on the adiabatic PES for 
oxygen dissociation on the usual simple metals such as Al(lll) and Mg(OOOl), obviously, our 
present prediction of the activated PES in the O2/Be(0001) system is unique. In particular, 
since beryllium has the same crystal structure and valence electrons as magnesium does, it 
is attempting to assume that the behavior of oxygen dissociation at Be(OOOl) is the same 
as or similar to that at Mg(OOOl). The result in Fig. 7, however, shows that this intuitive 
expectation needs to be reshaped. 

In the case of Al(lll) or Mg(OOOl), the failure of adiabatic DFT calculations in producing 
an activated-type PES has been ascribed to the unphysical output that charge transfer occurs 
even at large molecule-metal distance, which has led to speculations that nonadiabatic effects 



may play an important ro 



the simple sp bands 
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dissociation process at these metal surfaces with 
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23j. The unphysical charge transfer 



at large molecule-metal distance is generally caused by the requirement in the adiabatic 
description that the electron chemical potential of the O2 molecule aligns with that of the 
metal surface in the combined system (metal plus molecule). On the other hand, however, 
this requirement does not necessarily results in artificial charge transfer between the two 
largely-separated subsystems. Therefore, there is no universal criterion to judge whether 
the adiabaticity during molecular dissociation process breaks down or not. The molecule- 
metal interaction is so species-sensitive that even for the two metals with the same crystal 
structure and valence electrons, the molecular dissociation process at the surfaces can display 
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qualitatively different behaviors. To make it more sense, we have carried out a comparative 
study on the long-distance charge transfer effect by choosing the metal surfaces as Be(OOOl), 
Mg(OOOl), and Al(lll). The distances of the O2 molecule from the three kinds of metal 
surfaces are fixed at the same value of 7 A, at which charge transfer should not occur in 
reality. Also, the O2 molecules in the three systems are identically chosen to be in the 
HH-Z entrance. The calculated spin-split PDOS for the O2 molecule and the charge-density 
difference Ap(r) in the three systems are shown in Fig. 8. Clearly, one can see that for 
02/Al(lll) and O2/Mg(0001), although the distance of the O2 molecule from the surfaces 
is set to be as large as 7 A and charge transfer should not be presented in real systems, the 
adiabatic DFT calculations give yet the contrary results. In these two cases, it reveals in 
Fig. 8(b) and 8(c) that the LUMO of the ideal O2 molecule (i.e., the spin-down antibonding 
71* MO) has shifted down to align with the Fermi energy and becomes partially occupied, 
accompanying with an observable decreasing of the molecular spin. The 0-0 bond length 
also undergoes a little expansion compared to the free value of 1.235 A. Actually, Fig. 
8(b) and 8(c) depict nothing more than what have been known in previous studies, and 
it is exactly due to this unphysical (metal-molecule or intramolecule) charge transfer that 
motivates theoretical workers to remedy the calculations by, for example, nonadiabatically 
co„*ai„.ng the molecula. spins. As a d.ect .esult, a sizeable ene.gy ba.ie. w,U appear on 
the corrected PES [21]. When the attention is paid to the O2/Be(0001) system, however, 
we find that such unphysical large-distance charge transfer effect does not happen in the 
calculation. In fact, from Fig. 8(a) one can see that at a metal-molecule distance of 7 A, the 
calculated spin-split molecular PDOS shows no change with respect to the free case. The 
LUMO keeps empty and the charge- density difference is zero everywhere. Subsequently, 
the molecular bond length and spin of O2 are not influenced at all by the presence of the 
Be(OOOl) surface at a distance of 7 A. Combining with Fig. 7, therefore, we conclude that 
unlike Al(lll) and Mg(OOOl), the adiabatic DFT calculation presented in this paper is 
sufficient to predict an activated-type dissociation process of O2 at Be(OOOl). 

For further illustration of O2 dissociation at Be (0001), we plot in Fig. 9 four snapshots 
for the spin density (the contour spacing is 0.2 yU^/A'^) evolving along the dissociation path 
shown in Fig. 7. The values of the corresponding heights h and relaxed bond lengths d 
of the O2 molecule are also indicated in the figure. One can see from Fig. 9 that as the 
O2 molecule approaches the surface along the T-Y channel, the 0-0 bond length increases. 
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FIG. 8: (Color online) spin-polarized DOS for the molecular orbits of O2 with a fixed vertical 
distance of 7 A from (a) Be(OOOl), (b) Mg(OOOl), and (c) Al(lll) surfaces. The 0-0 bond length 
in each case is allowed to relax. The corresponding charge-density differences Ap(r) are plotted on 
the right side with the contour spacing of 0.002/A^. Solid and dotted lines denote the accumulated 
and depleted densities, respectively. 

while the spin magnetic moment decreases and tends to vanish at /i~2.3 A. Since there 
has occurred electron transfer from the substrate to the molecule at this height, thus the 
disappearance of the spin moment is due to the donation of the minority spins from the 
substrate not due to the spin flip inside the molecule. 

The adiabatic PES's for O2 dissociation at Be(OOOl) along the other parallel channels 
have also been calculated, which are found to have very similar elbow shapes except for 
visible differences in the transition states and dissociation barriers. To be more clear, a 
detailed comparison of these PES's is depicted in Fig. 10 by separately plotting the one- 
dimensional cuts of PES's and 0-0 bond length as functions of the O2 height h from the 
substrate surface. The insets in each panel in Fig. 10 display the initial and final positions 
of the two O atoms. Clearly, the dissociative adsorption of O2 along each of these parallel 
channels is a direct and activated type. The obtained dissociative adsorption energies, the 
dissociative energy barriers, and the geometrical parameters for the corresponding transition 
states, including the height of the O2 molecule and the 0-0 bond length, are concluded in 
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(a) h=3.9A 
d=1.236 A 




(b) h=3.1 A 
d=1.240 A 



(C) h=2.3A 
d=1.282 A 



(d) h=1.5A 
d=1.615 A 




FIG. 9: Snapshots of the spin quenching process during O2 dissociation at Be(OOOl) along the T-Y 
channel. The spacing of the spin-density contours is 0.2 ^b/^- The corresponding O2 molecular 
bond length d and distance h from the surface are also indicated. 
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FIG. 10: One-dimensional cuts of the potential-energy surfaces and the corresponding 0-0 bond 
lengths as functions of the O2 distance h from the surface, for eight different dissociative channels. 
The inset in each panel indicates the initial (filled circles) and final (hollow circles) atomic positions 
of O2. 
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TABLE II: The calculated dissociative adsorption energy {Ead), dissociative energy barrier (AE), 
and geometric parameters (namely, the O2 molecular bond length d and distance h from the surface) 
of the transition state, along the eight different dissociation channels. 



Channel Bad AE h (TS) do-o (TS) 



T-X 


4.00 


0.28 


2.10 


1.30 


T-Y 


4.14 


0.23 


2.10 


1.30 


B-X 


3.99 


0.26 


2.20 


1.29 


B-Y 


4.10 


0.26 


2.20 


1.29 


HH-X 


3.81 


0.35 


2.20 


1.30 


HH-Y 


4.09 


0.41 


2.10 


1.30 


FH-X 


4.00 


0.25 


2.20 


1.29 


FH-Y 


4.10 


0.25 


2.20 


1.28 



Table II. It is found that the dissociation path with both the lowest energy barrier and 
the largest dissociative adsorption energy is along the T-Y channel. For all the dissociation 
paths, the 0-0 bond lengths at the transition states are similar, ranging from 1.28 to 1.30 A, 
elongated from the 1.236 A in an isolated O2 molecule. Note that no dissociative adsorptions 
are found for the O2 molecules along the vertical entrances, which thus are not plotted in 
Fig. 9 and listed in Table II. This is different from the 02/Al(lll) system, in which the 
molecular dissociation along the vertical channels has been predicted to occur as well |l2|. 



VI. CONCLUSION 



In summary, by performing ab initio simulations we have for the first time systemat- 
ically investigated the adsorption and dissociation of oxygen molecules on the Be(OOOl) 
surface. We have identified both the physisorbed and the chemisorbed molecular states, 
which directly occur along the parallel and vertical channels, respectively. For the most sta- 
ble chemisorbed molecular state (which is along the HH-Z entrance), in particular, we have 
studied its electronic and magnetic properties by calculating the charge-density difference, 
the spin density, and the PDOS, which clearly show the charge transfer from the spin-up 
71* MO to the substrate followed by the back donation from the substrate to the spin-down 
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TT* MO, as well as the distinct covalent weight in the chemical bonding between O2 and 
Be(OOOl). This important covalent weight in the molecule- metal bond has been revealed 
through the combining fact: (i) The charge is largely accumulated along the Be-0 bond 
[Fig. 5(a)]; (ii) The a, c*, and tt MO's are heavily hybridized with the Be sp states (Fig. 
6); (iii) The chemisorbed molecular state in the HH-Z channel is much stable than that in 
the FH-Z channel, implying an observable role played by the subsurface Be atom beneath 
the O2 molecule. 

The energy path for the dissociation of O2 on Be(OOOl) surface has been determined 
by calculating the adiabatic PES's along various channels, among which the T-Y channel 
has been found to be the most stable and favorable for the dissociative adsorption of O2. 
Remarkably, our results have shown that the adiabatic dissociation process in the present 
O2/Be(0001) system is an activated type, with the lowest energy barrier of 0.23 eV in the 
most stable T-Y channel. To our knowledge, this is the first time to theoretically predict a 
sizeable adiabatic energy barrier during dissociation of O2 at the metal surfaces with simple 
sp bands. Thus as a final concluding remark, here we point out that in spite of many previous 
revealing and specific studies, an in-depth and primary insight into the common nature of 
the O2 dissociation at various simple sp metal surfaces remains yet to be attainable. 
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